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Abstract 

In this study we present two indicators that will reflect the difference be- 
tween collapsars and normal collapse-driven supernovae. They are products 
of explosive nucleosynthesis and neutrino emission. In the collapsar model, it 
is natural to consider that the product of explosive nucleosynthesis depends 
on the zenith angle because the system becomes highly asymmetric in order 
to generate a fire ball. We also consider the detectability of the HNRs which 
is located nearby our Galaxy. As a result, the number of the HNRs is esti- 
mated to be 5 X (10^ - 10~^), whose chemical composition can be spatially 
resolved. Using the optimistic estimate, more HNRs will be found and it will 
be possible to discuss on the chemical composition statistically. As for the 
energy spectrums of neutrinos, they are not thermalized in a collapsar be- 
cause the density of the accretion disk is much lower than that of a neutron 
star. The energy spectrums of (anti-)electron neutrinos from hypernovae will 
be mainly determined by the process of electron (positron) capture on free 
proton (neutron). It is also noted that high energy tail is not dumped in the 
case of hypernovae because the density of emitting region is low. Total energy 
of neutrino from hypernovae will depend on a lot of physical parameters such 
as total accreting mass and mass accretion rate, which are quite contrary to 
the situation of the normal collapse-driven supernovae. Therefore there will 
be a large variety of total neutrino's energies among collapsars. In the case of 
SN 1998bw, we think that the matter around the equatorial plane might be 
ejected from the system, which resulted in the formation of relatively weak 
jets and faint GRB 980425. 

PACS number(s): 98.70.Rz, 97.60.Bw, 97.10.Cv, UTAP-364, RESCEU-13/00 
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Typeset using REVT^ 



I. INTRODUCTION 



It will be a big progress that the fact that at least a part of the gamma-ray bursts (GRBs) 
comes from the hypernova (HN) explosions is being supported by the observations. It was 
reported for the first time that there seems to be a physical connection between GRB 980425 
and SN 1998bw [||]. They discovered an optical transient within the BeppoSAX Wide Field 
Camera error box of GRB 980425. Then they reported that the optical transient can be 
interpreted to be the light curve of SN 1998bw. As for the explosion energy of SN 1998bw, 
it was estimated to be as high as (20-50) x 10^^ erg as long as we believe the explosion is 
spherically symmetric 0,0]. This is the reason why SN 1998bw is called as a HN. The late 
afterglow of GRB 970228 also suggests the physical connection between GRB and HN 0. 
It was shown that the optical light curve and spectrum of the late afterglow of GRB 970228 
are well reproduced by those of SN 1998bw transformed to the redshift of GRB 970228. The 
afterglow of GRB 980326 is also believed to be the evidence for the GRB/HN connection 
due to the same reason 0. 

If we believe that a part of the GRBs comes from the explosion of massive stars, the 
explosion must be a jet-induced one because spherical explosion model has a difficulty in 
avoiding the baryon contamination problem ||^. In fact, some observations on GRBs are 
interpreted as evidence for the jet-induced explosion. For example, the breaks in the rate of 
decline of several afterglows can be explained by the beaming effect . The light curve and 
spectrum of SN 1998bw also seem to suggest a jet-induced explosion There are also 

some excellent numerical simulations on the jet-induced explosion of massive stars whose 
aim is to reproduce the fire ball [|10H T^, although the fire ball has not been reproduced yet. 

Here we must note the following two points, (i) it is not determined that all of the GRBs 
come from HNe. (ii) the explosion energy of SN 1998bw may be small if the explosion is 
the jet-induced one. Taking these points into consideration, we can classify the relation of 
GRB, SN, and HN as shown in Figure |l|. For example, region (a)/(d) means that HNe which 
didn't/ did generate GRBs. We note that HN fl SN = by definition. Here we defined that 
SN is the explosion of a massive star whose total explosion energy is about 10^^ erg. HN is 
defined as the explosion of a massive star whose total explosion energy is significantly larger 
than 10^^ erg. As for the region (b), other systems such as the merging neutron stars |[T3[ 
may belong to this region. 

One of the most famous model to realize a GRB from a death of a massive star is the 



collapsar model [|l4| , p!Tip!5| . The definition of the collapsar is written in |jl5[ as a massive 
star whose iron core has collapsed to a black hole that is continuing to accrete at a very 
high rate. Woosley also pointed out that there will be two types for collapsars. One (type 
I collapsar) is that the central core immediately forms a black hole with an accretion disk. 
The other (type II collapsar) is that the central core forms a neutron star at first, but the 
neutron star collapses to be a black hole with an accretion disk due to the continuous fall 
back. In both types, a strong jet, which is required to produce a GRB, is generated around 
the polar region due to the pair-annihilation of neutrinos that come from the accretion disk 
and/or MHD processes. The remnants of a collapsar will belong to the regions (a), (b), (c), 
and (d) in Figure 0. When the explosion energy of a collapsar is small, it will be classified 
as SNR. When the hydrogen envelope exists, a collapsar can not produce a GRB. 

Here we note that there are no observations that support directly the scenario of col- 



2 



lapsars. This situation is a contrast to that of the scenario of collapse-driven SN, which is 
supported by the detection of neutrinos at Kamiokande [16| and 1MB [17|. Thus we present 
in this study two observable indicators that reflect the mechanism of collapsars. These ob- 
servations will affirm the difference between collapsars and collapse-driven SN clearly. These 
are products of explosive nucleosynthesis and neutrino emission. We will discuss these es- 
sential features in the following sections. We also discuss the possibility of detection of such 
observations taking the event rate into consideration. 

In section |I|, we consider the explosive nucleosynthesis in the collapsar model. The 
luminosity and spectrum of neutrino from collapsars are shown in section |ITT[ Summary and 
discussion are presented in section |V|. 



II. FEATURES OF REMNANTS OF COLLAPSARS 



If we believe the collapsar model, the system is highly asymmetric in order to generate a 
fire ball PHHI^. Thus it is natural to consider that the product of explosive nucleosynthesis 
depends on the zenith angle (see details below). So we consider its detectabihty in this 
section. 

Here we must note that the asymmetric explosion can occur in the collapse-driven super- 
nova as long as the effects of rotation and/or magnetic field are taken into consideration, e.g. 
see p!8|^0|. As a result, the products of explosive nucleosynthesis also depends on the zenith 
angle in such an asymmetric collapse-driven supernova [0,^ . This means that we can not 
distinguish well whether it is the remnant of a collapsar or of a rotating collapse-driven 
supernova when we find an asymmetric SNR. In order to avoid such a problem, we should 
search for the hypernova remnants (HNRs), whose explosion energy can not be attained by 
the scenario of delayed explosion for the collapse-driven supernova. Thus we consider the 
detectabihty of the HNRs in this section. 

To tell the truth, the products of explosive nucleosynthesis in collapsars are not known 
exactly. There are many possibilities. For example, it is pointed out that ^^Ni is synthesized 
by the wind blowing off the accretion disk in a type I collapsar as long as the disk viscosity 
is set to be high [|l^ . In their simulations, the outflow containing much of ^^Ni is shown to 
be moving at 15 to 40 degrees off axis. However, the region where most of ^^Ni is contained 
may be around the polar region in a type I collapsar. This is the result of the explosive 
nucleosynthesis behind the strong jet. This picture is like the situation which occurs in 
the jet-like explosion of collapse-driven supernova ||2^. On the other hand, the chemical 
composition of the remnant of a type II collapsar may be spherically symmetric, because 
the launch of the jet is too late to cause the explosive nucleosynthesis . 



Such a situation is a contrast to that about the explosive nucleosynthesis in SN. The 
results of numerical calculations on explosive nucleosynthesis in collapse- driven SNe are 
compared with the observations very carefully, e.g. see ||2^p5[| . Thus, observations of the 
HNRs are necessary in order to determine which model is realistic and which model is 
unrealistic. Such observations may also give a light on the occurrence frequency of the type 
I collapsar relative to the type II collapsar. 

Here we estimate the chance probability to flnd the nearby HNRs whose chemical com- 
position can be resolved by the latest X-ray telescopes such as the ESA's X-ray Multi-Mirror 
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(XMM) and Chandra (AXAF) satellites whose spatial resolution is of order of 1 arcsec. 

At first, we consider the event rate of HN in a Galaxy. If we consider the event rate 
of GRB is equal to that of HN, the estimated HN rate becomes (10~^ - 10~®) yr~^ per 
Galaxy p6|-|28[|. If we take the beaming effect into account, the HN rate becomes larger 
than the observed GRB rate. On the other hand, the HN rate can be estimated to be 
~ 10^^ yr^^ per Galaxy, when we assume that the slope of the initial mass function is 
-1.35 the maximum mass of a star is 50Mq |^D|, (1O-3O)M0 stars explodes as collapse- 



driven SNe p5|], (3O-5O)M0 stars explodes as HNe, and the collapse-driven SN event rate is 



10~^ yr~^ per Galaxy This will be an upper limit for the HN event rate because all of 
the massive stars in the range (30-50) are assumed to explode as HNe. So we consider 
that the HN rate is in the range (10~^ - 10~^) yr^^ per Galaxy. 

In order to know the chemical composition of the ejecta, HNR must be so young that the 
remnant is not composed mainly by the inter-stellar medium (ISM) but by the HN ejecta. 
Here we consider the Fe distribution in the remnant because the main products of explosive 
nucleosynthesis, ^^Ni, decays to Fe. We can estimate the shock radius, Rg, at which the 
amount of Fe from the ejecta becomes equal to that from ISM in the remnant as follows: 

where Mpe, n, /i are the amount of Fe from the ejecta, mean ambient hydrogen density, and 



mean atomic weight of cosmic material per H atom [32|. Here we assumed that the mass 



fraction of Fe in the ISM is equal to that in the solar system abundances I^S]. On the other 



hand, the shock radius in the adiabatic phase can be written as follows 




Rs = 7.9 ^ —- [pc], (2) 



where i^exp and t are total explosion energy and age of the remnant, respectively. In the 
case of SN 1998bw, Mpe is estimated to be 0.7 Mq 0. If we consider that this is the standard 
case with HN, t has to be less than the following value: 

That is, roughly speaking, the HNR whose age is less than 10^ yr must be searched for in 
order to know the chemical composition of the ejecta. As for the limit of the distance from 
the Earth to the target, it must be nearer than 3 Mpc in order to resolve the asymmetry of 
the chemical composition of the ejecta as long as the spatial resolution of the X-ray telescope 
is 1 arcsec. 

Since there are 55 galaxies within 3 Mpc from our Galaxy [Q, the number of the HNRs 



is estimated to be 5 x (10^ - 10"'^), whose chemical composition can be spatially resolved. 
(10~^ - 10~®) yr~^ per Galaxy and 10^ yr are adopted for the HN event rate and the age of 
the oldest HNR, respectively. Using the optimistic estimate, the HNRs will be found more 
and we will be able to discuss on the chemical composition statistically. 



Here we consider the report of Wang on NGC 5471B and MF83 in MlOl. They 



reported that NGC 5471B and MF83 may be the HNRs since they require explosion energies 
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comparable to the energies frequently associated with GRBs. Since the distance of MlOl 
from our Galaxy is about 7.2 ± 0.4 Mpc |^|, it seems difficult to observe the distribution 
of the chemical composition of the remnants at a present state. However, we can say that 
the HNR event rate seems larger than the lower estimate for the ORB rate if these are 
really HNRs. Although other interpretations are possible for these highly luminous X-ray 
sources ||3^ , search for the hypernova remnants nearby our Galaxy has a potential to reveal 
the mechanism of the GRB. 



III. NEUTRINO EMISSION FROM COLLAPSARS 



The second important feature of coUapsars is that no neutron star but an accretion 
disk around the black hole is formed. One of the most probable heating source for the jet 
formation is believed to be the vV annihilation emitted from the accretion disk [^. On the 
other hand, neutrinos are emitted from only the surface of a neutron star in SN explosion. 
In this section we discuss the differences of the energy spectrum of the emitted neutrinos 
between the collapse-driven SN and the collapsar. 

In the case of SN, the energy spectrum of neutrinos is approximately represented by the 
thermal distribution, because the mean free path of neutrinos is much shorter than the radius 
of the neutron star . Strictly compared to the perfect Fermi-Dirac distribution with zero 



chemical potential, however, the high energy phase space is less populated [40,|T|. This is 
because the high energy tail is dumped due to much larger opacities (oc el). In addition, 
it is well-known that the total energy of neutrinos is determined by only the gravitational 
binding energy of the neutron star p9|| . 

On the other hand, the energy spectrums of neutrinos are not dumped in coUapsars, 
because the nucleon density of the accretion disk is much lower than that of a neutron 
star [Q. Namely the energy spectrums of neutrinos emitted from coUapsars are entirely 
proportional to the emission rates. Then the total energy of neutrinos could depend on a 
lot of physical parameters such as the total accreting mass M, the mass accretion rate M, 
and so on. Therefore there will be a variety of total energies of the emitted neutrino for 



coUapsars. MacFadyen and Woosley |11| have shown that the accretion disk in a collapsar 



can be described well by the analytic solution derived by Popham et al. WM. According to 



their analytic solution, the neutrinos are mainly emitted from the region where T = (1-10) 
MeV and p = (10^-10^°) g cm'^. 

If we assume that the density and the temperature are constant in the neutrino emitting 
region , we can estimate the energy spectrum of the emitted neutrinos from the accretion 
disk. For n + ^ p + i^, the spectrum of z7g in unit time, unit volume, and unit energy is 
represented by 



j2 eAf r<2 



where E^^ is energy of z7g, T is temperature, Gp is Fermi coupling constant, Ca — 1-37 is 
normalized by the experimental value of neutron lifetime Tn — 887.6 s [^], is number 
density of neutron, Q ~ 1.29 MeV, and rrie is electron mass. For e~ + z/g + 14? "we 

obtain 
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^2 e+e- ^2 i .qo ( ^2 _ ( /T^2\3/2 



where Cy = 1/2 + 2sin^ Ow, Ca = 1/2, and sin^ 9w — 0.231 is Weinberg angle and we 



assume E^^ ^ T. 

In Figure |^(a) we plot the obtained spectrum of z7g emitted from the accretion disk 
{d^riy^l dtdEi^^ = d^rfj^ / dtdEy^ + d^iifjj^^ /dtdE^J in unit time, unit volume, and unit energy. 
It should be noted that the high energy tail is not dumped at all because the nucleon density 
of the accretion disk is much lower than that of a neutron star and the mean free path is 
much longer. This is remarkable feature only for coUapsars. 

The luminosity of z7g can be obtained by integrating the spectrum as g = 
/ dE^^Ey^d'^n^J dtdEiy-^. Then we obtain q"^ ^ 4.6 x lO^^pioTf^Xnucerg cm"^ s"^ and 
qe+e- ^ 2.4 X lO^^Tf^erg cm'^ s'^, where pio = p/10^° g cm'^, Tn = T/10^^ K, and X^^^ 
is the mass fraction of nucleons. Xnuc is given by Xnuc = 30. 97p7o^''^Ti^o^exp(— 0.6096/ Tio), 



where X„nc < 1 44 



As is clear from the above relations, the luminosity of neutrinos from collapsars de- 
pends sensitively on the temperature. Therefore if the configuration of the accretion disk 
is modified by the change of the environment such as the mass accretion rate, mass of the 
progenitor, and the mass of the black hole, then the total luminosity and energy of neutrino 
from collapsars will be changed drastically. These points are entirely different from the nor- 
mal collapse-driven SN because the total energy of neutrinos from SN is determined only by 
the gravitational binding energy of the central neutron star. The event numbers of z7g from 
HN expected at Super-Kamiokande is represented by 

dR VaNp d^n^-^ rj^ \a+ fR\ 

Vp^-(Eg+)-— — (Ee+)At, (6) 



dE^+ AttD^ 'dtdE^-^ 

where Ef.+ = Ey^ — Q is the energy of the positron which is scattered through p + Ve n + e~^ 
in the detector, apij^ = 2^(1 + 3(7^ )Ef,+ J El+ — ml is the cross section of the process, Va is 
the volume of the emitting region in accretion disk, A'p ~ 1.5 x 10'^^ is the number of proton 
in Super-Kamiokande, D is the distance from the earth to the coUapsar, At ~ M/ M is the 
duration of the emission. In Figure |^(b), we show the plot of the event number adopting an 
representative parameter set We can find that the neutrino emission from a collapsar 
can be observed at Super-Kamiokande as long as it is located within ~3 Mpc from the earth. 
As for the estimate of the detection rate at Super-Kamiokande becomes as follows: 

P ~ 5 X (10^^ - 10"^) [yr-^], (7) 

where the same way of estimation is done as in section That is, (10~^ - 10^^) yr^^ 
per Galaxy and 55 are adopted for the HN event rate and the number of galaxies within 3 
Mpc from our Galaxy. Using the optimistic event rate ~ 5 x 10~^ per year, the detection 
probability of the collapsar can be as large as that of the collapse-driven SN. 



IV. SUMMARY AND DISCUSSION 
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In this study, characteristic products of nucleosynthesis and neutrino emission have been 
proposed as two indicators that will reflect the features of the collapsars. 

We consider the detectability of the HNRs because we can not distinguish well whether 
it is the remnant of a collapsar or of a rotating collapse-driven supernova when we find 
an asymmetric SNR. As a result, the number of the HNRs is estimated to be 5 x (10^ - 
10"^), whose chemical composition can be spatially resolved. Using the optimistic estimate, 
more HNRs will be found and it will be possible to discuss on the chemical composition 
statistically. Due to such observations, we will be able to determine which model is realistic 
and which model is unrealistic. Such observations may also give a light on the occurrence 
frequency of the type I collapsar relative to the type II collapsar. Moreover, we can say that 
the HNR event rate seems larger than the lower estimate for the GRB rate if NGC 5471B 
and MF83 in MlOl are really HNRs. Although other interpretations are possible for these 
highly luminous X-ray sources search for the hypernova remnants nearby our Galaxy 



has a potential to reveal the mechanism of the GRB. 

Strictly speaking, there will be a little difference between the SNRs of collapsars and 
those of collapse-driven supernovae. We think that an extreme jet-induced explosion like 
collapsars will not happen in the case of SN. This is because almost all of the matter has 
to be ejected in order not to leave a black hole but to leave a neutron star at the center. 
That is, matter around the equatorial plane has to be also ejected, which will be observed as 
'jet-like' explosion like SN 1987A ||22|. On the other hand, an extreme jet-induced explosion 
is required in order to make fire balls for the model of the jet-induced HN. So, even if the 
matter around the equatorial plane is ejected due to some reasons in the case of HN too, 
the degree of jet-induced explosion will be very large and chemical composition will depend 
strongly on the zenith angle in the case of the type I collapsar. 

It is also noted that the mass accretion rate becomes low if the matter around the 
equatorial plane is ejected from collapsars. This will result in the decline of the total energy 
of neutrinos emitted from the accretion disk. As a result, total explosion energy may become 
small in that case. It is reported that the explosion energy of GRB 980425, which is said 
to be associated with SN 1998bw, is quite lower than that of the usual GRBs [Q. In the 
case of SN 1998bw, we think that the matter around the equatorial plane might be ejected 
from the system, which resulted in the formation of relatively weak jets and faint GRB 
980425. This means that SN 1998bw and GRB 980425 may be classified in the region (e) in 
Figure [l|. Of course, this picture requires that the system of SN 1998bw and GRB 980425 
is highly asymmetric, because the total explosion energy of SN 1998bw is estimated to be 
(20-50) X 10^^ when spherical explosion is assumed 

As for the (anti-)electron neutrino emission from the collapsars, its energy spectrum 
is mainly determined by the emission rate due to electron (positron) capture on proton 
(neutron). As the temperature becomes higher, contribution of the process of electron- 
positron pair annihilation can not be negligible. It is also noted that high energy tail is not 
dumped in the case of the collapsar because the density of emitting region is low. These 
features on energy spectrum are quite different from that of SN. 

Total energy of neutrino depends on many physical quantum such as total accreting 
mass and mass accretion rate. It is noted the emission rate due to the electron capture on 
proton is proportional to T^. So a little change in temperature results in great change in 
the neutrino flux. That is why there will be a variety of total luminosity of neutrino among 
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collapsars, which is in striking contrast to the case of SN. As for the event rate, the detection 
probabihty of the collapsar can be as large as that of the collapse-driven SN if we use the 
optimistic event rate ~ 5 x 10~^ per year at Super-Kamiokande. 

Finally, we stress again that these features on nucleosynthesis and neutrinos will reveal 
the mechanism of GRB quite well. We hope the increase of further observations in the near 
future. 
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FIGURES 




FIG. 1. Classification of the relation of GRB, SN, and HN. Here we defined that SN is the 

explosion of a massive star whose total explosion energy is about 10^"*^ erg. HN is defined as the 
explosion of a massive star whose total explosion energy is significantly larger than 10^"*^ erg. As 
for the region (b), other systems such as the merging neutron stars (Ruffert & Janka 1999) may 
belong to this region. 



11 




FIG. 2. (a) Energy spectrums of iJ^. from collapsars. (b) Event numbers expected at Su- 
per-Kamiokande. Solid line represents the total energy spectrum. Dashed line represents the 
contribution from n + e^ ^ p + z7g. Dotted line represents the contribution from e"'" + e~ — > i/g + z^e- 
The temperature, the nuclear density, and the volume of the emitting region are set to be 5 MeV, 
]^q10 g (.j]^-3^ g^j^^ 6. 5x10^*^ cm'^, respectively (Popham et al. 1999). The distance, the total accret- 
ing mass, and the mass accretion rate are set to be 3 Mpc, SOMq, and O.IMq s~^, respectively. 
Then the total event number which is obtained by dE^-\- integration is ~ 15. 
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